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SIGNAL COMPRESSION BASED ON CURVATURE PARAMETERS 



Related Applications 

This document is related to co-pending and commonly assigned U.S. Patent 
Application Serial No. 09/703,269, entitled "CURVATURE BASED METHOD 
FOR SELECTING FEATURES FROM AN ELECTROPHYSIOLOGICAL 
SIGNALS FOR PURPOSE OF COMPLEX IDENTIFICATION AND 
CLASSIFICATION," inventor Sweeney et al., filed on October 3 1 , 2000, the 
specification of which is hereby incorporated by reference in its entirety. 

This document is related to co-pending and commonly assigned U.S. Patent 
Application Serial No. 10/195,838, entitled "USE OF CURVATURE BASED 
FEATURES FOR BEAT DETECTION," inventor Sweeney, filed on July 15, 2002, 
the specification of which is hereby incorporated by reference in its entirety. 

This document is related to co-pending and commonly assigned U.$. Patent 

e 

Application Serial No. XX/YYY,YYY, entitled TACHYARRHYTHMIA 
DETECTION AND DISCRIMINATION BASED ON CURVATURE 
PARAMETERS," inventor Sweeney et al., filed on June 27, 2003, the specification 
of which is hereby incorporated by reference in its entirety. 

Technical Field 

This document relates generally to signal compression and particularly, but 
not by way of limitation, to electrocardiogram signal storage in an implantable 
medical device. 

Background 

Implantable cardiac rhythm management devices typically monitor and 
process cardiac signals to provide therapy to the heart. Therapy may include 
delivering a pacing pulse to trigger a contraction of the heart or delivering a 
defibrillation shock to interrupt an abnormal heart rhythm. Some cardiac rhythm 

Attorney Docket No. 279.635US1 ] 



management devices also monitor cardiac performance or other physiological 
parameters for use in controlling the delivery of pacing pulses. 

There exists an unmet need for providing improved data storage systems, 
devices, and methods to provide better diagnosis or treatment of patients. 

Summary 

'i 

The present subject matter provides a curvature based method of selecting 

-w 

features from a sampled signal. In one embodiment, the sampled signal includes a 
cardiac signal or an electrocardiogram. According to one aspect of the present 
subject matter, characteristic points are selected by sensing a cardiac signal on a 
real-time basis and upon receipt of each sample, computing curvatures on a 
continuous basis. Characteristic points are sometimes referred to as significant 
points or data points. 

In one embodiment, each "turn" in the signal is represented with a 
characteristic point that has measures of time, value, area and width. Time refers to 
a time denoting the center of the turn in the curvature series. Value refers to the 
amplitude of the input signal at the time of the center of the turn. Area refers to a 
value denoting the direction and degree of the turn. Width refers to a value denoting 
the duration over which the turn occurs. 

In one embodiment, the area is described as the angle taken by the signal 
between the start of a turn and the end of a turn where the turn defines a lobe. The 
area, in one embodiment, is expressed in radians. 

By the methods presented herein, each new characteristic point is known as 

soon as the turn which it describes has ended. In one embodiment, the stream of 

..... ... ... 3 < 

characteristic points are analyzed continuously as each new characteristic point 
occurs. In one embodiment, the characteristic points are saved into a buffer for 
analysis at a later time. 

All turns in an input signal are not the same. Some turns may be large 
deflections with sharp angles (for example, turns in a QRS signal) while other turns 
might be slight and not associated with big deflections (for example, noise in the 
signal). In one embodiment, the present subject matter includes a scheme to select 
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prominent characteristic points in a signal and discard the smaller ones that may 
represent noise or only slight signal deflections. 

The area for the characteristic point is an indirect measure of the angular 
turn in the signal. The maximum turn is approximately 180 degrees which 
represents a very rapid signal reversal. A large deflection in the signal, such as the 
R-wave (a full signal swing in 10 samples over, a period of 50 ms) would have large 
turn angle of 178 degrees. A small noise spike (such as a 10% full scale swing for a 
single sample) may also have a large turn angle of 175 degrees. The present subject 
matter treats these signal deflections differently. The R-wave has a large area 
whereas the noise spike would have a small area since it occurred over a lower 
number of samples. Thus, the characteristic point area is a surrogate for 
determining which characteristic points are important and which are not. 

The effect of the foregoing is that small characteristic points at or near the 
signal baseline value are ignored. According to one embodiment, an absolute area 
threshold value in the range of 0.5 is used for a rate-sensing channel. 

In one embodiment, the present subject matter automatically adjusts the 
threshold value to better adapt to varying patients and signals. 

In one embodiment, each characteristic point includes a time component, an 
amplitude component and an area component. According to the present subject 
matter, the time, amplitude and area components are derived from a curvature 
calculation performed using the sampled signal. 

In one embodiment, the plurality of characteristic points are stored in lieu of 
the sampled signal. By storing characteristic points rather than the sampled signal, 
the memory required for a particular duration of a sampled signal is reduced, or 
compressed. Further compression of the data is achieved, in one embodiment, by 
storing a repetition marker denoting that a group, or cluster, of characteristic points 
is substantially a duplicate of another cluster of characteristic points. Each cluster, 
for example, can represent a particular morphological feature such as a heart beat. 
A timing code is stored with the repetition marker to denote when the cluster 
occurred. In one embodiment, a difference marker is stored to indicate that a cluster 
of characteristic points differs from another cluster of characteristic points. A 
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difference code is stored with the difference marker to denote differences between 
the particular cluster and a previously stored cluster of characteristic points. As 
with the repetition marker, a timing code is stored with the difference marker to 
denote when the particular cluster occurred. 
5 In one embodiment, only characteristic points having an absolute value of an 

area larger than a predetermined threshold are stored. Characteristic points having 
an absolute value of the area less than the threshold value are discarded or ignored. 

In one embodiment, the sampled signal is reconstructed from the stored 
characteristic points by establishing a function between sequential characteristic 
10 points. In various embodiments, the function includes a straight line, a polynomial 
fit, a parabolic fit, a cubic spline or other function. 

Other aspects of the invention will be apparent on reading the following 
detailed description of the invention and viewing the drawings that form a part 
thereof. 

15 

Brief Description of the Drawings 

In the drawings, like numerals describe substantially similar components 
throughout the several views. Like numerals having different letter suffixes 
represent different instances of substantially similar components. 
20 Fig. 1 illustrates an arbitrary signal as a function of time. 

Fig. 2 illustrates osculating circles fitted to an arbitrary signal as a function 
of time. 

Fig. 3 illustrates a signal expressed in analog-to-digital converter value as a 
function of sample number. 
25 Fig. 4 illustrates a signal expressed in analog-to-digital converter value as a 

function of sample number and having a constant curvature. 

Fig. 5 illustrates characteristic points based on curvature as a function of 

time. 

Fig. 6 graphically illustrates a method for determining area under a curvature 

30 lobe. 
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Fig. 7 illustrates tabulated data for compressing as a function of a repeated 
cluster of characteristic points. 

Fig. 8 illustrates tabulated data for compressing as a function of a cluster of 
characteristic points that differs from a previous cluster of characteristic points. 
5 Fig. 9 illustrates a method to determine a characteristic point according to 

one embodiment of the present subject matter. 

Fig. 1 0 illustrates a method to find and store repeated groups of 
characteristic points according to one embodiment of the present subject matter. 

Fig. 1 1 illustrates a block diagram of an implantable device according to one 
10 embodiment of the present subject matter. 

Detailed Description 

In the following detailed description, reference is made to the accompanying 
drawings that form a part hereof, and in which is shown, by way of illustration, 

1 5 specific embodiments in which the invention may be practiced. These embodiments 
are described in sufficient detail to enable those skilled in the art to practice the 
invention, and it is to be understood that the embodiments may be combined, or that 
other embodiments may be utilized and that structural, logical and electrical changes 
may be made without departing from the scope of the present subject matter. The 

20 following detailed description is, therefore, not to be taken in a limiting sense, and 
the scope of the present invention is defined by the appended claims and their 
equivalents. 

This document discloses, among other things, systems, method and 

■ 

apparatus involving cardiac rhythm management systems used for sensing a cardiac 
25 signal, including intracardiac electrogram or surface electrocardiogram (ECG). 
Such cardiac rhythm management systems include, but are not limited to, 
pacemakers, CRT devices, cardioverter/defibrillators, pacer/defibrillators, and drug 
delivery devices. However, it is to be understood that the present methods and 
apparatuses of compressing a signal can be applied to heart beat detection as well as 
30 other signals related to cardiac activities, including, but not being limited to, 
mechanical motion, sound, pressure, acceleration, and impedance signals. 
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Curvature-Based Analysis 

For an arbitrary voltage versus time signal, the salient characteristics can be 
captured in terms of those points along the signal where it makes characteristic 
5 "turns." In general, these points are different from those that would be selected 
using first or second derivative criteria as shown in Fig. 1 . In the figure, which 
illustrates arbitrary electrogram 100, the R wave, denoted as point 101, has a high 
second derivative d V/dt . Further, the electrogram sections on either side of point 
101 have high positive and negative slopes detectable with a first derivative criteria. 
10 However, those points with these large (or even maximum) slopes do not convey 
any particularly significant description of the electrogram. For example, at each 
time along a segment between points 101 and 102, the waveform has a large 
negative slope but no point along this segment stands out significantly from any 
other. 

15 On the other hand, points 102, 103, 104 and 105 are neither maximums nor 

minimums of the electrogram or its derivatives. These points are descriptive of the 
arbitrary waveform shape. Points 102, 103, 104 and 105 are salient because they 
mark locations where the signal makes significant turns. The turn at point 101 is 
very sharp and the turns at points 103 and 105 are less sharp turns and more broad. 

20 The turns at points 102 and 104 are even less sharp but rather local. The present 
subject matter detects points 101, 102, 103, 104 and 105 using a criteria based on 
the curvature of the signal. 

According to one embodiment of the present subject matter, signal curvature 
is illustrated in Fig. 2. In the figure, arbitrary electrogram 120 is shown. At each 

25 point along electrogram 120, an osculating circle can be found that fits tangentially 
with the local portion of the electrogram. Each circle has a radius. The curvature of 
electrogram 120 at that point is the inverse of that radius so that small circles have 
large curvatures while large circles have small curvature. Fig. 2 shows these circles 
at selected points. At point 122, the circle is larger, and thus the curvature is 

30 smaller, than at points 121, 123 and 124. At points 121, 123 and 124, the turns are 
sharper and the curvature is larger. 
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In general, the curvature at point (X, Y) of an arbitrary curve in a two- 
dimensional space is expressed as: 

Curvature = d 2 Y/dX 2 

[l + {dY/dX} 2 ] 3/2 

5 

As seen, the curvature is a non-linear combination of both the first and 
second derivatives of the curve. At those points along the curve where the first 
derivative is zero (for example, point 101 in Fig. 1), the curvature is equal to the 
second derivative and at points where the second derivative is zero (for example, 
10 any straight section regardless of its slope), the curvature is zero. 

The present subject matter calculates curvature of an arbitrary input signal 
on a sample-by-sample basis. 

Consider the question of dimensionality of the curvature for electrogram 
signals. When both X and Y have the same dimensions (say length) then curvature 
15 has a dimension of 1 /length. For an arbitrary signal having a voltage as a function 
of time V(t), such as an electrocardiogram, the signal is transformed into a time- 
versus-time signal T(t) according to T(t) = V(t)/U where U is a constant with 
dimensions of voltage/time. With this transformation, the first and second 
derivatives of T(t) become 



20 



dT_(t) = dV_(t) • 1 / U which is dimensionless and 

dt dt 



25 d 2 T (t) ~ d V (t) « 1/U which has dimensions of 1/time 



dt 2 dt 2 



and thus the curvature has dimensions of 1/time. Curvature is then expressed as: 

30 Curvature = d 2 T(tVdt 2 = d 2 V(tVdt 2 /U 

[ 1 + {dT(t)/dt} 2 ] 3 1 2 [ 1 + {d V(t)/dt / U} 2 ] 3 ' 2 

which has dimensions of 1/time and U has a numerical value. 

Consider next, curvature as a function of signal gain or amplitude of the 
35 input signal. Assume arbitrary gain G is applied to the input signal to find a new 
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input signal F(t) wherein F(t) ~ G «T(t) = V(t) • G/U. The curvature of the gained 
signal is then: 

Curvature = d 2 F(tVdt 2 = d ■ VftVdt 2 » G/U . 

[1 + {dF(t)/dt } 2 ] 3/2 . [1 + {dV(t)/dt . G/U} 2 ] 3/2 

5 

Since the selection of the gain applies to all point in the input signal, the gain 
G can be thought of as a value that emphasizes the curvature at different points 
along the signal. For a given V(t), a small value for G will emphasize the curvature 
at those points on V(t) where dV(t)/dt is large and will deemphasize the curvature at 

10 those points where dV(t)/dt is small. Large values for G have the opposite effect. 

The ratio G/U can be expressed as W having dimensions of time/voltage. In 
one embodiment, the input signal is a voltage sampled with a 1 2-bit analog-to- 
digital converter (ADC) having numerical voltage values in the range +/- 2048 
where each value represents a number of basic amplitude units AV, or voltage 

1 5 resolution. It is assumed that the amplifiers in the present subject matter are 
adjusted so that the samples from V(t) largely fill this range. Furthermore, it is 
assumed that the samples are taken at a fixed rate and thus the time is represented by 
an integer number of samples with each sample representing a time interval AT = 
time resolution = 1 /(sample rate). 

20 Fig. 3 illustrates numerical X-Y space 130 which is denoted by sample 

numbers on the X-axis and the analog-to-digital conversion value on the Y-axis. 
Consider square 135 in X-Y space 130 which is 256 steps along the X-axis and 256 
steps along the Y-axis and rectangle 140 in voltage-time space 145 that this square 
represents. Rectangle 140 has a width of 256/(sample rate), in seconds and a height 

25 of 256/(voltage resolution), in volts. Depending on the values selected for AV and 
AT, this rectangle may or may not be a square in voltage-time space. 

In one embodiment, W is selected to require that a square in voltage-time 
space be represented by a square in sample-sample space. Under that condition, the 
curvature versus time relationships that exist in voltage -time space are preserved in 

30 sample-sample space. 

In Fig. 4, an artificial voltage-time curve is shown that forms a semicircle 
inside square box 150. Traversing from left to right across box 150 in voltage time 
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space, the curvature of the signal is zero until the semicircle is encountered. The 
curvature jumps to a constant value, equal to 1 /radius of the circle, and again jumps 
to zero at the end of the semicircle. For a particular value of W, the representation 
of this signal in sample-sample space 155 also has a constant curvature and for other 
5 values of W, the curvature would not be constant in sample-sample space. 

To maintain this relationship between voltage-time and sample-sample 
space, W is selected as follows. In voltage-time space, the box T wide by V high is 
taken to be square. In sample-sample space, the box is T/AT time-samples wide and 
GV/AV voltage-samples high and through a transformation, the voltage sample is 

10 converted into time-samples using U so that the sample-space square has 

dimensions of T/AT by VG/AV/U. For the box to be square in sample-space, and 
assuming the box (V by T) in voltage-time space is square, then (WAV) • G/U =. 
T7AT or G/U = W - (T/AT)/(V/AV). 

According to one embodiment of the present subject matter, curvature is 

15 based on the first and second derivative of the signal. A least square cubic 

polynomial fit is used to reduce the noise that would otherwise result from using 
numerical estimates for the derivatives and using non-linear calculations to find 
curvature. 

The sampled voltage signal is expressed as V(t) - V(I ♦ AT) where tH • AT 
20 so that the fit of size N uses 2N+1 voltage samples centered on the time t, thus: 

■ 

V([i-n].AT), . . .V([i-2].AT), V([i-1].AT), V(i.AT), V([i+1].AT), V([i+2].AT), . . 
.V([i+n].AT). 

25 This set of N sampled data points is used to make the least-squares cubic fit 

which is given as V est [i.AT + dt] - Ai + Bi • [dt] + Ci • [dt] 2 + Di . [dt] 3 where Ai, 
Bi, Ci, and Di are coefficients determined by minimizing the square error for the fit 
and dt represents a time step away from i»AT at which V es t is evaluated. The 
coefficients in the polynomial are denoted with an i to indicate that they are valid 

30 for the point i AT. 

Using the above equations, the curvature on a point to point basis, at the time 
UAT, becomes 
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Curvature ( i.AT) - d 2 V(tVdt 2 . W = 2 v Ci . W . 

[1 + {dV(t)/dt.W} 2 ] 3/2 [1 + {Bi.W} 2 ] 3/2 

Sample points, however, will not necessarily fall at the times where the 
signal curvature has a maximum or minimum value. Thus, in one embodiment, the 
5 curvature signal is integrated between adjacent sample points using estimates for the 
first and second derivative of the signal at the sample points. 

Further simplification yields an expression for average point curvature as: 

'. Cavgi = 1 • IC^B, 2 -OB^ B;,, -C ; +0 , + OBu 2 -Oi+B^Bi +C; -Cl . 1 , Bj B M 

(Bj-Bj i) 2 7 (B 2 +l) ,/2 (B M U l) r/2 " 

10 

Cav gi ^ ± • Ci + Ci-i , Bj =B,! 
2 (Bi 2 + 1) 3/2 

As noted, curvature is computed on a sample-by-sample basis from the input 

15 signal. 

Consider next, a procedure for finding the characteristic points in the signal. 
Turns in the original signal are reflected as excursions above and below zero 
in the curvature signal. As suggested by curvature 160 of Fig. 5, each lobe above 
zero (for example, lobe 165) or below zero (for example, lobe 170) then represents a 

20 single turn in the input signal. Positive curvature lobes represent a leftward turn and 
negative curvature lobes represent a right ward turn. The area under each lobe is the 
total angle included in the turn. 

A point-by-point method is used to identify the lobes as they occur and to 
find the area and centroid of each lobe. 

25 Curvature values are generated on a sample by sample basis. At each 

sample time, CRV represents the current curvature value and CRV 0 | d represents the 
previous value that is retained from the previous sample. When the curvature value 
is zero, the original signal is not turning so that no characteristic points may exist. 
The value CRV will rarely be exactly equal to zero. 

30 In one embodiment, a dead-zone surrounding zero is defined where the 

calculated curvature may be treated as equal to zero. For times where the curvature 
is within the zone, the signal is not changing significantly. The zone is defined with 
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a single curvature threshold value that extends above and below zero from 

+CRV t hresh tQ -CRVthresh. 

In one embodiment, the present subject matter identifies nine cases for 
considering the value of CRV relative to the threshold and the absence or direction 
of a lobe. These cases can be described as follows: 

Case 1 : CRV>CRV t hresh and not in a lobe. 

Here, the current curvature value is above the dead-zone and the curvature 
signal is not currently in a lobe. Thus, a positive lobe has just started so the positive 
lobe initialization calculations will follow. 

Case 2: CRV t hrcsh >CRV > -CRVthresh and not in a lobe. 

Here, the current curvature values is inside the dead-zone and the curvature 
signal is not currently in a lobe. 

Case 3 : CRV<- CRVthresh and not in a lobe. r 

Here, the current curvature value is below the dead-zone and the curvature 
signal is not currently in a lobe. Thus, a negative lobe has just started so the 
negative lobe initialization calculations will follow. 

Case 4: CRV>CRV t hresh and in a positive lobe. 

Here, the current curvature value is above the dead-zone and the curvature 
signal is in a positive lobe. Thus, the positive lobe continuation calculations will 
follow. 

Case 5: CRVthresh >CRV>- CRVthresh and in a positive lobe. 

Here, the current curvature value is inside the dead-zone and the curvature 
signal is in a positive lobe. Thus, a positive lobe has just ended so the positive lobe 
finalization calculations will follow. 

Case 6: CRV< -CRVthresh and in a positive lobe. 
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Here, the current curvature value is below the deadrzone and the curvature 
signal is currently in a positive lobe. Thus, a positive lobe has just ended and a 
negative lobe has just started so that both the positive lobe finalization calculations 
and the negative lobe initialization calculations will follow. 

Case 7: CRV>CRV t hresh and in a negative lobe. 

Here, the current curvature value is above the dead-zone and the curvature 
signal is in a negative lobe. Thus, a negative lobe has just ended and a positive lobe 
has just started so that both the negative lobe finalization calculations and the 

positive lobe initialization calculations will follow. 

/- 

Case 8: CRV t hresh>CRV> -CRVthresh and in a negative lobe. 

Here, the current curvature value is inside the dead-zone and the curvature 
signal is in a negative lobe. Thus, a lobe has just ended so the negative lobe 
finalization calculations will follow. 

Case 9: CRV< -CRVthresh and in a negative lobe. 

Here, the current curvature value is below the dead-zone and the curvature 
signal is in a negative lobe. Thus, the negative lobe continuation calculations will 
follow. 

In one embodiment, hysteresis is used to aid in identifying curvature lobes. 
After having started a lobe, one embodiment provides that the curvature value cross 
a threshold value closer to zero in order for the lobe to finish. Thus, hysteresis 
introduces another threshold value. 

Consider next, the metrics for characterizing each lobe according to the 
present subject matter. In one embodiment, those metrics include the total area of 
the lobe, the time of the lobe area centroid and the value of the original data at the 
time of the area centroid. In one embodiment, other metrics or parameters are used 
to describe the characteristic point, including, for example, the peak curvature in the 
lobe, the time of the peak curvature, the times of the lobe start or lobe finish, a time 
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of the centroid, a peak, a midpoint, a midpoint at a median value or any other 
predetermined parameter. 

Fig. 6 shows an example of curvature lobe 185 as a series of calculated 
curvatures at each sample time. Tick marks 1 80 along the X axis represent the 
actual samples and thus are separated by AT. Zero curvature is shown as the X axis. 
The threshold curvature values, -^CRVthresh 190 and -CRVthresh 195 are shown as 
horizontal lines above and below the X-axis. Curvature points 200-245 are shown 
as small solid circles. 

As shown in the figure, initial curvature values at curvature points 200 and 
205 are within the dead-zone of the threshold curvature values and thus, no lobe is 
yet established. Curvature point 210 is above the threshold and corresponds to Case 
1 described above. The magnitude of initial area 250 is calculated as soon as the 
lobe is started according to case 1 described above. Upon determining curvature 
point 215, a contribution to the magnitude of main area 255 is calculated according 
to case 4 described above. In addition, the magnitude of main area 255 is increased 
with the determination of curvature points 220, 225, and 230. Upon determining 
curvature point 235 lying below +CRV thresh 190, the magnitude of final area 260 is 
calculated for the region below curvature lobe 185 and between curvature point 230 
and the intersection with +CRV t h re sh 190. Curvature points 235, 240 and 245 lie in 
the dead-zone and do not contribute to an area calculation. The area of curvature 
lobe 185 is calculated to include the sum of initial area 250, main area 255 and final 
area 260. 

In one embodiment, a value M of a counter is maintained by a processor to 
monitor sample-by-sample development of a lobe. Initially, it is assumed, the signal 
will be in the dead-zone and for a sample excursion outside the dead-zone, the value 
M is reset to zero and incremented with each subsequent point that remains outside 
the dead-zone. Values for M are illustrated at the different curvature points in Fig. 
6. 

To find the area under the curvature lobe, the value of M is set to zero and 
the initial area of a curvature lobe is computed as initial area = (Ci + CRVthresh) • (Ci 
- CRVthresh) / (Ci-Cj.i) when the lobe is started. For each successive curvature value, 
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the value of M is incremented and the additional area contribution is computed as 
additional area - - (M-l) • Q + M • Cm. When the lobe ends, the value of M is not 
incremented and the final area contribution is computed as final area = M • Ci-i + 
(Cj + CRVthresh) # (Ci - CRVthresh)/(Ci-Ci-i). The sum of the initial, final, and all 
5 additional areas may then be multiplied by AT/2 to find the area under the curvature 

m 

lobe. 

To find the first moment of the area under the curvature lobe, a similar 
approach is used. When the initial area of the curvature lobe is computed, the initial 
moment of the area is ajso computed as initial moment = - {C\ + 2 CRVthresh} • y 2 

10 where y - (CRVthresh - CM)/(Ci-Q-i). For each successive curvature value in the 

lobe, an additional moment contribution is computed using the same M values as for 
the area computations. This additional moment contribution is computed as 
additional moment = {Cj • (3M-1) + Ci-i • (3M-2)}. When the lobe ends, a final 
moment contribution is computed at the same time that the final area is computed. 

15 The final moment is computed as final moment -[CRVthresh • { 3M + 2y} + Cm • 
{3M + y}]*y. 

The sum of the initial, final, and all additional moments may then be 
multiplied by AT 16 to find the first moment of the area under the curvature lobe. 
The time of the centroid of the curvature lobe area is found by dividing the first 
20 moment of the curvature lobe area by the curvature lobe area as time ~ first moment 
of area/area. This time is with respect to the time of the curvature point that started 
the lobe. 

For both the area and the first moment of the area, the calculations are 
constructed in such a way that the contribution of each curvature sample may be 
25 summed into a total area or total first moment of the area as the samples are 

collected. Thus, for all calculations, the factors include the present and previous 
curvature values, the threshold values, and the counter M. In this manner the area 
and the time of the centroid of the area are generated as soon as the lobe is ended. 

In one embodiment, the present subject matter identifies characteristic points 
30 generated from point-by point processing of an input signal. In one embodiment, 
the characteristic points in the signal are detected and saved into a buffer. In one 
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embodiment, each characteristic point is a set of values including the time of the 
characteristic point, the value of the input signal at that time, and a value describing 
the direction and extent of curve in the signal that produced the characteristic point. 
In one embodiment, each time a heart beat is detected, the characteristic 
5 points for that beat are removed from the buffer and compared to characteristic 
points of one or more previously stored beat types, thus classifying each beat. 

The beat classification process determines if the set of characteristic points 
describing the present unknown beat matches with the set of characteristic points 
describing a previously known beat type. 

10 In one embodiment, if a match is found, then a repetition code is stored 

indicating that the current set of characteristic points is a repeat of a previously 
stored beat. The repetition code includes a value for the time of the occurrence of 
the current beat and includes identification of the previously stored beat. In one 
embodiment, the repetition code identifies the previously stored beat and a code 

15 denoting the number of repeated beats and a number corresponding to the time 
interval between each successive beat. 

In one embodiment, a difference code is stored to indicate that the current set 
of characteristic points differs from the set of characteristic points of a previously 
stored beat. The difference code includes a value for the time of the occurrence of 

20 the current beat and includes identification of the previously stored beat as well as 
the specific difference(s) from the previously stored beat. 

Consider next a method of extracting characteristic points. In one 
embodiment, characteristic points of an electrocardiogram are extracted using a 
sampling at 400 Hz. In one embodiment, sampling is performed at 200 Hz. The 

25 sampled data is filtered using a 5-point or 3-point running average filter and a 5 
least squares regression average-point curvature method to select characteristic 
points relating to the QRS complex. 

For detecting curvatures associated with slower morphologies, sampling is at 
a slower rate, such as, for example, 50 Hz. An example of a wave having slower 

30 morphologies would be the T-wave in an electrocardiogram. 
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In one embodiment, multiple sampling rates are used for characteristic point 
extraction. 

In one embodiment, dual-rate sampling is performed with sampling at 
200Hz providing fast characteristic points and sampling at 50 Hz providing slow 
5 characteristic points and a 3 -point running average filter is used. 

In one embodiment, curvature is found by using 5-point least squares 
regression filtering. The coefficients for finding the linear (Bi) and parabolic (Ci) 
fit-coefficients in the cubic least squares regression to the data at point i are as 
follows: 

10 





Linear 




Parabolic 


P-2= 


(1/12)* Rate 


Q-2 


= (10/70)* Rate 2 


P-l- 


(-8/12) * Rate 


Q-! 


= (-5/70)* Rate 2 


Po = 


0 


•Qo 


= (-10/70) * Rate 




(8/12) * Rate 


Qi 


= (-5/70) * Rate 2 


P 2 = 


(-1/12) * Rate 


Q 2 


= (10/70)* Rate 2 



15 

where Rate is the sampling rate in sample/sec, Bi is the sum of the products of these 

P least squares regression coefficients multiplied by the 5 corresponding data points 

centered on i (Di_2 to Dj +2 ) and Ci is the sum using the Q least squares regression 

20 coefficients. Curvature is then computed as: 

Curvature = 2 * Ci * W 

[1 + {Bi.W} 2 ] 3/2 

where W is a constant. Since W appears in conjunction with fit-coefficients Ci or 

25 Bi, the value can be incorporated into the computations as a data signal gain such 

that D'(i) = W • D(i). Then, the curvature expression becomes 

Curvature = 2 • C [ i 

[l + IB'i} 2 ] 372 

30 where B'i and Ci are found using the D' data points rather than the D data points 
with the least squares regression coefficients. 

The expression for curvature then becomes: 

CurV = Curvature . Gn = C ' i . 2 . Gn . = C"i 



[1 + {B'i} 2 ] 3/2 [i+ {B'i} 2 ] 372 



35 
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where Gn is an arbitrary value which is incorporated (with 2) into C M . If the 
curvature thresholds used for detection of characteristic points is adjusted 
accordingly, it does not matter whether point-by-point curvature or Gn times the 
point-by-point curvature is calculated. Incorporating the values 2 • Gn into C"i can 
5 be accomplished by changing the least squares regression coefficients Q_ 2 through 
Q2 into new least squares regression coefficients (Q'_2 through Q2 ) by 
multiplication with 2Gn. If Gn is selected as 7/(6 Rate) and the term Fn = Rate/12 is 
incorporated as additional data signal gain, then the least squares regression 
coefficients become: 



10 Linear Parabolic 

P"- 2 = 1 Q M -2= 4 

P"_i = -8 Q"_i= -2 

P"o =0 Q" 0 = -4 

P" , = 8 Q M , ~ -2 

15 P" 2 = -1 Q M 2 = 4 



with D"(i) = 2 . W . Fn ♦ D( i ); and with Ci" = Q". 2 • D"(i-2) + Q"-i • D"(i-1) + Q M o • 
D'*(i) + Q M + i . D M (i+l) + Q\ 2 • P"(i+2); and with Bi" - P n _ 2 • D M (i-2) + P M _i . 
D n (i-l) + PVD M (i) + P"-i •D ,, (i+1) + P" +2 • D M (i+2). 

Note that the coefficients are all powers of 2 (either 0, 1, 2, 4 or 8), thus 
20 simplifying hardware or firmware implementation of the present subject matter. 

When using use two sample rates for determining characteristic points, the 
computed curvatures will be different when computed at different rates. Thus, two 
threshold levels and parallel sets of computations are used. Also, the running 
average filtering is different with the fast rate filtering using 5 consecutive points 
25 while the slow rate filtering uses every fourth point spaced over 16 points. 

In one embodiment, the linearity of the least squares regression operation is 
used. The running average filtering operation are performed at the fit-coefficient 
level (i.e. the Bs and Cs) rather than the data level. In this embodiment, two 
filtering operations are performed (one for Bs, one for Cs) at each rate. 
30 In one embodiment, the fast and the slow running average filtering 

operations are performed in parallel. 

To classify the beats for compression purposes, the following 
implementation is used to determine the characteristic points. 
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Starting at the time of one sample, the sensing hardware and software of the 
present subject matter begins acquiring the next sample so that the actual sample 
reflects the nature of the signal between the last sample and the current sample. At 
the actual time of the sample, the value of the sample may be ascribed to the time 
5 half-way between the current sample arid the previous sample. Samples at the slow 
sample rate apply to the signal two fast sample steps earlier that the actual sample 
time. 

In one embodiment, a circular data buffer is used. Data is acquired at the 
fast sample rate and used to fill a circular data buffer. Once the timing of a fast or 

10 slow characteristic point is determined, the data value for that characteristic point is 
found from the data stored in the circular buffer. Thus, the data buffer is sized so 
that the requisite data is still available in the buffer. In one embodiment, the buffer 
is sized to hold about 500 millisecond of data or 256 data points. In one 
embodiment, the circular buffer is implemented as a memory array with an index 

1 5 pointer. In one embodiment, the index pointer includes a binary counter and the 
buffer size is a power of 2. 

In one embodiment, running average filtering (RAF) is performed on the 
input data. In various embodiments, the filtering is provided by a circular buffer or 
a shift register set. The values in the circular data buffer reflect the fast running 

20 average filtering. In one embodiment, a hybrid approach is used where the fast 
running average filtering is performed and then the data is placed into a circular 
buffer. In one embodiment, the running average filter includes a shift register 
approach. The fast sample running average filtering registers are called DFreg o to 
DFreg -3. In one embodiment, an accumulator SUM p is also maintained. 

25 When a new fast sample (D 0) is obtained, the following operations are 

carried out: 

SUM F =SUM Fold + Do 
SUM Fold =SUM F - DFreg. 1 
DFreg .3 = DFreg .2 
30 DFreg -2- DFreg _i 

DFreg -1 = DFreg 0 
DFreg 0 = D 0 
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The value of SUM F is the sum of the 3 most recent data points. Thus, 
circular buffer CircBuf is filled. Next the CircBuf [CBIndex] is set to SUM F? and 
CBIndex is incremented followed by resetting CBIndex to zero on overflow. The 
value SUM p is also loaded into the fast data register set for the least squares 
regression operation. 

On every 4 sample at the fast rate, a slow rate sample is obtained as 
SLOWVAL= SUM F + DFreg -3 and the slow running average filtering operation is 
performed using the following: 

SUM S = SUM sold + SLOWVAL 
SUMsoid = SUM S - PSreg., 
DSreg .1 = DSreg 0 
DSrego^ SLOWVAL 

With the above operations, the circular buffer is loaded with the results of 
the fast running average filtering as the sum of the last 3 samples. The average 
values are determined by dividing by 5 (or 3). The fast running average filtering 
average applies to the signal at a time 1 .5 time steps back from the time of the 
current sample. 

On each 4 m sample, SLOWVAL is the sum of the previous 4 fast samples 
and thus applies to the signal 2 time steps back from the current sample. The result 
of the slow running average filtering operation is the sum of the last 3 SLOWVAL 
samples and thus tbe sum of the last 12 fast samples. Again the average is found by 
dividing the sum by 12. The slow running average filtering average applies to the 
signal 6 time steps back from the time current sample. 

To perform the least square regression fit (cubic fit) operations, the filtered 
data point is multiplied by W • Fn where Fn is equal to the sample rate divided by 
12. In one embodiment, the slow rate is selected to be 1/4 of the fast rate (RATE) 
and the fast rate and slow rate samples are multiplied by W • Fnf = W ♦ RATE/ 1 2 and 
W • Fn s = W • RATE/48 = W • Fnf I A respectively. In one embodiment, the value of 
W is selected such that the product W • Fn f is a power of 2 which also makes W ♦ Fn s 
a power of 2. 
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In one embodiment, the data used for least squares regression computation is 
extracted from the circular buffer by manipulating pointers into the buffer to select 
the correct values. In one embodiment, the data is shifted through a set of 5 
registers and the fit computations are performed on the register set. 

The fast set of shift registers are called Freg o through Freg ,4 and the slow 
registers are called Srego through Sreg .4. The shift operations are: 

Freg .4 = Freg .3 Sreg .4 ~ Sreg _ 3 

Freg _ 3 - Freg _ 2 . Sreg -3 = Sreg _ 2 

Freg _ 2 Freg _ { Sreg _ 2 - Sreg _i 

Freg -i - Freg 0 Sreg -1- Srego 

and the data is loaded into the register sets as Freg 0 = SUMf • W • Fnf and Sreg 0 = 
SUM s .W.Fn,/4. 

In one embodiment, the fast register set load/shift/computation operations 
occur on each sample while the slow register set load/shift/computation operations 
occur on every 4 sample. 

In one embodiment, for fast-sampling fit, the least squares regression fit 
coefficients are obtained by B"f = (Freg .4 - Freg 0) + 8 • (Freg _i - Freg .3); and C" F 
= 4 • (Freg .4 + Freg 0 ) - 2 • (Freg _i + Freg .3) - 4 • Freg -2. The linear and parabolic 
fit coefficients for the fast sample rate are denoted B M F and C M F respectively. 

At 200 Hz sampling, the time in the signal at which these cubic least squares 
regression fit coefficients apply is 3.5 time steps earlier than the time of the current 
sample. At 200 Hz, 3-point running average filtering is used and the sample again 

■ 

reflects the signal 1/2 a time step back from the time of the current sample. 
Furthermore, the running average delay is 1 step, the least square regression delay is 
2 steps. 

For slow sampling fit, the linear and parabolic fit-coefficients for the slow 
sample rate are denoted B M s and C"s respectively and are found from the slow 
register set as B M s = (Sreg .4 - Sreg 0) + 8 • (Sreg_i - Sreg-3); and C'V = 4 • (Sreg-4 + 
Srego) - 2 • (Sreg .1 + Sreg. 3 ) - 4 . Sreg. 2 . 

The mechanics of filtering, regression, and curvature calculation are the 
same as the fast sample rate. The time in the signal at which these slow cubic least 
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squares coefficient fits apply is 3 1/2 slow (14 fast) time steps earlier than the time 
of the current sample. Each slow sample is the sum (or average) of the 4 recent fast 
samples. Thus, when a slow sample occurs (i.e. on each 4 th fast sample), the slow 
sample has the same time as the concurrent fast sample time. However, the slow 
sample applies to the signal 1/2 slow time steps (i.e. 2 fast time steps) earlier than 
the current slow sample time. The 3 -point running average delay is 1 slow (4 fast) 
samples and the least squares regression delay is 2 slow (8 fast) samples for a total 
of 3 1/2 slow (14 fast) sample time steps. 

In one embodiment, a fast curvature value is computed for each current fast 
sample from the fast linear (Bf!) and parabolic (Cf M ) least square regression fit 
coefficients. In one embodiment, a slow curvature value is computed for each 
current slow sample from the slow linear (B s ") and parabolic (C S M ) least square 
regression fit coefficients. 

In one embodiment these curvatures are computed as point curvatures as 
Curvature = C7 (1 + B") 3/2 where C" and B M are the current fast or slow fit 
coefficient. 

In one embodiment, these curvatures are computed as average-point 

* 

curvatures as 

Curvature 1 (C n "B" 2 -C"B"B o " ^C" -KV + C"B o " 2 -C n ,, Bn M B"+C"-C n M l; B" f B n " 

(B"-B 0 ") 2 " (B" 2 +l) ,/2 . " (B 0 2 + 1) 1/2 

Cavgi ~- J_ • 9l + C o" , B M = B 0 M 
2 (B ,,2 + l) 3/2 

where C" and B" are the current fast or slow fit coefficients and Co" and Bo" are the 
previous value for these same coefficients. 

In one embodiment, running average filtering reduces curvature noise 
relative to the signal content of fast curvature signals. An electrocardiogram signal 
typically takes relatively large excursions away from zero curvature at the times of 
fast turns in the signal. However, when the fast turns are less severe, the curvature 
signal does not move as far away from zero and curvature noise can prematurely 
end the curve lobe. The result is that one small curvature lobe is broken into two or 
more even smaller lobes. Thus, in one embodiment, the curvature noise relative to 
the signal content is reduced, for fast curvature signals, by running average filtering. 
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Filtering of the input signal does not achieve the same effect as filtering the 
actual curvature signals. In one embodiment, 3 -point running average filtering is 
applied to the fast curvature values. 

In one embodiment, fast curvature values are computed as average point 
5 curvature and 3 point running average filtering of the curvatures so that the fast 
curvature value computed when the current sample is taken applies to a time in the 
input signal that is 5 time steps earlier than the current sample. 

In one embodiment, slow curvature values are computed at 1/4 the fast 
sample rate and are computed as average-point curvatures so that when the every 
1 0 fourth fast sample is taken, slow curvature is computed and it applies at a time in the 
input signal that is 4 slow (16 fast) time steps earlier than the current sample. 

In one embodiment, the detection of a curvature-lobe and characteristic 
points entails both fast and slow characteristic point detection performed in parallel. 
In one embodiment, separate curvature thresholds and curvature area thresholds are 
1 5 used for the fast and slow operations. 

In one embodiment, a curvature lobe is ended when the curvature signal 
crosses the same threshold used to start the lobe but in the opposite direction of the 
crossing that started the lobe. In one embodiment, a hysteresis is used so that a 
curvature lobe is ended when the curvature signal crosses a threshold that is smaller 
20 (i.e. closer to zero) than the threshold used to start the lobe but in the opposite 

direction of the crossing that started the lobe. In one embodiment, a curvature lobe 
is ended when the curvature signal crosses zero as the hysteresis value. 

In one embodiment, the value of W is selected according to the following 
procedure. Select an approximate range for W from the following analysis. At a 
25 200 Hz sampling rate, the QRS takes approximately 5 samples. The signal at the 
peak is approximated by 128. If the QRS were to be approximated (in curvature- 
space) as a semicircle arching from 0 to the peak then back to 0 in 5 points, then the 
radius of that semicircle would be 0.01 sec and the peak would have a value of 0.01 . 
Thus, the value of W would be the one that, when multiplied, the peak (128) would 
30 have the value of 0.01, i.e. W - 0.01/128 = 0.000078. In one embodiment, W = 
0.000122= 1/8192 which is a power of 2. 
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Consider selection of curvature threshold values for characteristic point 
detection. Thresholds, in one embodiment, reduce or eliminate the noise 
characteristic points and computational burden. This noise is the natural results of 
the curvature signal hovering around zero (or some other value) when the raw signal 
5 does not otherwise have characteristic turnings. 

In one embodiment, for a curvature lobe to be detected, the point- wise 
curvature falls outside of a dead-band formed by the threshold values on either side 
of zero. Reducing this threshold towards zero increases the number of curvature 
lobes that are detected and increases the computational burden. In one embodiment, 

10 a threshold value is selected by using approximately a 1 minute epoch of the input 
signal to make a histogram of the fast curvature values computed at each sample 
time. The threshold value is selected as one which corresponds to 20% of the peak. 
This value is then used as the threshold for both the fast curvature and slow 
curvature feature selection process. 

15 In one embodiment, the noise lobes are removed by requiring that a detected 

lobe must exceed a critical area limit. The noise lobes have generally small areas 
while the real lobes associated with characteristic turns in the signal usually have 
substantial areas. In one embodiment, the area limit of 0.1 is used. 

• ■ • ' , - r * ... ........ 

In one embodiment, different values are used for the fast curvature and slow 
20 curvature area thresholds. In one embodiment, these thresholds are pre-selected. In 
one embodiment, the thresholds are adapted to patient signals. Adapting the 
thresholds may occur once in a learning session or may occur periodically. In one 
embodiment, the thresholds are adapted by using the number of beats occurring in a 
one minute epoch to set the gain and curvature thresholds and selecting a desired 
25 number of fast and slow characteristic points to find for these beats. 

For example, in one embodiment, 5 fast characteristic points per beat and 8 
slow characteristic point per beat are used as targets. The area thresholds are then 
adjusted until these target numbers of characteristic points are found. In one 
embodiment, a lower area limit of 0.1 is used. 

30 
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Signal Compression 

In one embodiment of the present subject matter, data corresponding to the 
sampled signal is stored in the form of characteristic points. For example, rather 
than storing all data points from a sampled input signal, one embodiment provides 
5 that the characteristic points are stored. Each characteristic point, in one 

embodiment, includes a time of occurrence of that characteristic point and the 
amplitude of the input signal at that time. In one embodiment, each characteristic 
point includes data as a function of the size of the lobe at that time. 

In one embodiment, a minimum absolute value of a size of a characteristic 

10 point is established as a threshold value: In such a case, a characteristic point is 
discarded or not stored if the absolute value of the size is less than a minimum 
absolute value and stored if the value is greater than the minimum absolute value. 
For example, in one embodiment, a minimum size of 0. 1 is established. 

In one embodiment, the threshold value is stored in a memory accessible to a 

15 processor of the present subject matter. In one embodiment, the threshold value is a 
fixed value. In one embodiment, the threshold value is a variable and is adjusted 
manually or automatically. In one embodiment, the threshold value is varied as a 
function of the input signal or other measured parameter. For example, in one 
embodiment where the input signal is an electrocardiogram, the threshold value is 

20 varied as a function of the beat frequency. 

Other compression systems are also contemplated for use with repetitive 
signals. For example, in one embodiment where the input signal is an 
electrocardiogram, a library of beat types are classified and stored in memory. Each 
beat is represented by a series of characteristic points and using methods described 

25 herein, the set of characteristic points are assembled in a group or cluster. In one 
embodiment, subsequent beats, having a relationship with the previously stored 
group of characteristic points, are stored by referencing the previously stored beat 
group. 

In one embodiment, the relationship is one in which the subsequent beat 
30 matches a previously stored beat. Fig. 7 illustrates tabulated data for compressing 
as a function of a repeated cluster of characteristic points. In the figure, 
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electrocardiogram 295 includes beats 300, 310 and 320. In the figure, beat 300 is 
indicated to have occurred at time 0, beat 3 10 at time x\ and beat 320 at time x\ + 
Aj. Table 1 identifies cluster 305, cluster 315 and cluster 325, associated with beat 
300, beat 310 and beat 320, respectively. Each cluster includes a set of 
characteristic points labeled CP1-CP7 with each characteristic point having an 
associated time, amplitude and size. In memory, the characteristic points of cluster 
305 are associated with a particular beat type occurring at a specified time. 
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Table 1 

As shown in the figure, beat 310 and beat 320 are found to be a match with 
beat 300. Criteria for determining a suitable match between beats entails, among 
other things, establishing alignment between the beats and comparing relevant 
features^ 

In one embodiment, a repetition marker is stored in memory to indicate the 
occurrence of beat 310 and beat 320. The repetition marker, in one embodiment, 
includes identification of beat 300. The repetition marker, in one embodiment, 
includes timing information to denote that beat 3 1 0 occurred at time xj and that beat 
320 occurred at time Xi + Aj. In one embodiment, the repetition marker indicates 
that an interval of time between beat 310 and beat 320 is Ai . 

In one embodiment, a series of matched repeated beats is stored in memory 
as a duplicate of a previously stored beat along with a uniform interval of time 
between consecutive beats. 
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Fig, 8 illustrates exemplary electrocardiogram 345 having beat 350, 360 and 
370, each of which differ. Table 2 identifies cluster 355, cluster 365 and cluster 
375, associated with beat 350, beat 360 and beat 370, respectively. Each cluster 
includes a set of characteristic points labeled CP1-CP7 with each characteristic 
5 point having an associated time, amplitude and size. In this example, beat 350 

occurred at time t=0 and has a cluster of characteristic points CP1-CP7 as shown in 
cluster 355. Beat 360 occurred at time t= X2 and has a cluster of characteristic 
points as shown in cluster 365. Beat 370 occurred at time t- x 2 + A2 and has a 
cluster of characteristic points as shown in cluster 375. The entries in cluster 365 

10 for CP7 differ from those of CP7 of cluster 355. The entries in cluster 375 indicate 
that the size of each characteristic point is different from the size of each 
characteristic point in cluster 355. In one embodiment, beat 3 60 and beat 3 70 are 
stored in memory using difference markers. In one embodiment, the difference 
marker includes identification of beat 350. The difference marker, in one 

1 5 embodiment, includes timing information to denote that beat 360 occurred at time 
t= X2 and that beat 370 occurred at time X2 + A2. The difference marker, in one 
embodiment, indicates the differences in characteristic point CP7 of cluster 365 
relative to that of cluster 355. For example, according to one embodiment and with 
respect to the figure and Table 2, the difference marker for cluster 365 indicates 

20 CP7a = 51; 25; -58. The difference marker, in one embodiment, indicates the 
differences in size of each characteristic point of cluster 375 relative to that of 
cluster 355. For example, according to one embodiment and with respect to the 
figure and Table 2, the difference marker for cluster 375 indicates SIZE A = +10; 
-30; +40; 0; 0; -30; 90. Other encoding schemes are also contemplated for both the 

25 difference marker as well as the repetition marker. 
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Table 2 

Fig. 9 illustrates method 400 according to one embodiment of the present 
5 subject matter. At 405, a sample point is received. The sample point corresponds to 
data from a signal to be compressed. In one embodiment, a sample series is stored 
in memory and receiving the sample point entails accessing the memory. At 410, a 
value for the curvature is calculated for the sample point. Curvature is a function of 
the first and second derivatives at the sample point, as described elsewhere in this 

10 document. At 415, the series of curvatures values are processed to identify lobes 

above or below zero where each lobe corresponds to a single characteristic point. In 
one embodiment, the calculations to determine each characteristic point value is 
performed as the series of curvature values are processed as described elsewhere in 
this document. At 420, an inquiry is made to determine if a lobe is already started 

1 5 but not yet finished. If a lobe is not in progress, at 425 an inquiry is made to 

determine if the curvature sample is large enough to start a lobe and if so, at 430 the 
calculations for starting a lobe are made. If at 420 a lobe is in progress, at 435 an 
inquiry is made to determine if the curvature value is such that the current lobe has 
ended. If not, at 440 the calculations for the lobe are continued. If the lobe has 

20 ended, at 445 the calculations to end the lobe are made and at 450 the characteristic 
point values are found. At 455, an inquiry is made to determine if the characteristic 
point is large enough to be saved and if so, at 460 the characteristic point is saved 
either as the compressed signal or stored in a buffer for additional compression. It 

V 
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will be recognized that the above sequence of steps entails identifying a 
characteristic point. 

Fig. 10 illustrates method 500 for additional compression according to one 
embodiment of the present subject matter. At 505, the next characteristic point is 
5 received and placed into a buffer. At 5 1 0, an inquiry is made to determine if a beat 
has occurred in the signal to be compressed. In one embodiment this determination 
is made from the characteristic points in the buffer. In one embodiment this 
determination is made from other processing of the signal to be compressed. In one 
embodiment, the beat is a periodic component of the signal to be compressed. In 

10 one embodiment, the signal is a cardiac signal and the beat is a heartbeat. 

If a beat has occurred, at 515 one or more groups of characteristic points in 
the buffer are identified. In one embodiment, groups are identified based on the 
timing of characteristic points relative to the times of the current and/or previous 
beats. In one embodiment, the group includes all members in the buffer since the 

1 5 previous beat. In one embodiment, the groups are formed using the values of the 
characteristic points stored in the buffer. 

At 520, an inquiry is made for each identified group of characteristic points 
to determine if the group of characteristic points matches with a previously stored 
group. If the group does not match, then at 525 the group and time are stored in the 

20 compressed signal and processing returns to receiving the next characteristic point. 
If the group does match, then at 530, an inquiry is made to determine if the match is 
exact. If the match is exact, then at 535 an identifier for the matched group is stored 
in the compressed signal and at 540, the time for the group is stored in the 
compressed signal and processing returns to receiving the next characteristic point. 

25 If the match is not exact, then at 545, the differences from the matched group are 
found and at 550 the differences are stored and then at 535 and 540 the matched 
beat identifier and time are stored and processing returns to receiving the next 
characteristic point. 

In one embodiment, a library of previously stored groups is compared with 

30 each subsequent group to determine if a match exists. In one embodiment, a library 
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of previously stored groups is compared with each subsequent group to identify 
differences. In one embodiment, both differences and matches are identified. 

Reconstruction 

5 A reconstruction of the input signal can be generated as a function of the 

stored characteristic points. In one embodiment, for example, the input signal is 
modeled as a series of straight line segments connecting characteristic points in 
chronological order. In this embodiment, the area, or size, of the characteristic point 
is ignored or discarded and a plot is created using the time and amplitude, 

10 In one embodiment, a polynomial function is fitted to a series of 

characteristic points to reconstruct the input signal. In various embodiments, the 
polynomial function includes a parabolic function or a cubic spline function. In one 
embodiment, the area of the lobe is used in reconstructing the signal. In one 
embodiment, the start time for each lobe is used in reconstructing the signal. Other 

15 criteria can also be used for reconstructing the signal. For example, the area, the 
end time of the lobe, a peak value, or other parameter describing the lobe is stored 
and later used to reconstruct the signal. 

» 

Exemplary Device 

20 In one embodiment, the present subject matter includes circuitry, hardware 

and software for implantation in a body. In one embodiment, the present subject 
matter is adapted for use external to a body. 

Fig. 1 1 includes a diagram of one embodiment of portions of cardiac rhythm 
management system 600 according to the present subject matter. As illustrated, 

25 system 600 includes unit 603, lead 610 and lead 640, each of which are coupled to 
heart 605. In various embodiments, unit 603 includes an implantable unit and 
operates as a pacemaker, a CRT device, heart failure or cardiac resynchronization 
device, a cardioverter/defibrillator, a pacer/defibrillator and a drug delivery device. 
In the figure, unit 603 senses cardiac activities from heart 605 and delivers therapy 

30 to heart 605 according to programming and circuitry of unit 603. Leads 610 and 
640 provide electrical connections between unit 603 and heart 605 and each 
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includes an electrode adapted to be disposed in or about heart 605. Lead 610 is 
coupled to input circuit 615 and provides an electrical signal sensed from cardiac 
activity of heart 605. Lead 640 provides therapy to heart 605 as a function of a 
signal received from therapy circuit 635. 
5 In one embodiment, unit 603 includes input circuit 615, controller 620, 

telemetry circuit 625, memory 630 and therapy circuit 635. Input circuit 61 5 is 
coupled to lead 610 and receives an analog signal based on cardiac activity of heart 
605. In one embodiment, input circuit 615 includes an analog-to-digital converter. 

In one embodiment, controller 620 includes a processor and programming to 
10 implement a method as described herein. In one embodiment, controller 620 
includes a circuit to control therapy circuit 635 and telemetry circuit 625. 

Memory 630 is coupled to controller 620 and provides storage for 
compressed and uncompressed data. In one embodiment, the input signal is 
sampled and stored in memory 630 and at a later time, the data is accessed, 
15 compressed and stored again in memory 630. In one embodiment, programming 
executing on a processor of controller 620 is stored in memory 630. 

In one embodiment, controller 620 determines what therapy is to be 
delivered by therapy circuit 635 to heart 605 using lead 640. 

In one embodiment, controller 620 determines what data is to be 
20 communicated using telemetry circuit 625. Telemetry circuit 625, in various 
embodiments is adapted to provide wired or wireless telemetry. 

In one embodiment, controller 620 includes a comparator which generates 
an output based on a comparison of an input signal and a reference signal. 

L 

Exemplary comparisons include determining if the lobe is making a positive or 
25 negative excursion from a baseline. In one embodiment, an exemplary comparison 

includes determining if the curvature series lobe is greater or less than a dead-zone 

threshold. In one embodiment, an exemplary comparison includes comparing the 

curvature signal with a hysteresis value of a threshold. 

In the figure, two leads are shown coupled to heart 605. In one embodiment, 
30 a single lead provides both sensing and therapy for heart 605, 
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Alternative Embodiments 

Variations of the above embodiments are also contemplated. For example, 
in one embodiment, the present subject matter is adapted to compress data that is 
received from an implanted medical device. In one embodiment, the present subject 
5 matter is used to compress templates and waveforms for storage in an implanted 
medical device. 

In one embodiment, the present subject matter is used to compress an 
arbitrary waveform or signal. In addition to compressing voltage as a function of 
time, other parameters are also candidates for compression. For example, in various 
10 embodiments, a temperature signal, a sensor signal, an accelerometer, an impedance 
measurement or other varying signal is compressed using the present subject matter. 
In various embodiments, a digital or analog signal is compressed using the present 
subject matter. 

In one embodiment, a slowly changing signal, such as temperature as a 
15 function of time, QRS width as a function of time or QT as a function of time, is 
compressed using the present subject matter. 

The characteristic points may result from a fast sampling rate or a slow 
sampling rate. In one embodiment, a number of sampling rates greater than two are 
provided, wherein a sampling rate is stored with each characteristic point. In one 
20 embodiment, the threshold values are adjusted based on the sampling rate. 

In one embodiment, the time resolution for sampling is adjusted to achieve a 
desired compression ratio. In one embodiment, the signal amplitude resolution is 
adjusted to achieve a desired compression ratio. For example, with a slowly varying 
input signal, a lower resolution of the signal amplitude yields satisfactory results. In 
25 one embodiment, data denoting the time resolution and the signal amplitude 
resolution is stored with each characteristic point. 

In one embodiment, the present subject matter includes a data storage 
medium, such as a computer disk, with programming to execute a method as 
described herein. ' 
30 In one embodiment, the present subject matter is implemented in dedicated 

circuitry to perform a method described herein without use of a microprocessor and 
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programming. The circuitry may include hardwired analog or digital circuitry 
including, for example, operational amplifiers or logical gates. 

In one embodiment, the present subject matter stores a group of 
characteristic points associated with a particular pattern of excursions in the input 
signal and stores the group in memory. When the pattern occurs at a later time in 
the input signal, as identified by a subsequent group of characteristic points, a 
repetition, or similarity, marker is stored in memory. When a different pattern 
occurs at a later time in the input signal, as identified by a subsequent group of 
characteristic points, a difference marker is stored in memory. In one embodiment, 
the particular pattern includes one or more heart beats in an electrocardiogram. In 
one embodiment, the particular pattern includes, for example, a series of five regular 
heart beats followed by a spike. 

In one embodiment, the present subject matter includes a circuit to receive a 
sampled electrocardiogram signal and generate a plurality of characteristic points. 
As a function of the archival history of characteristic points stored in a memory, the 
present subject matter, in one embodiment, is programmed to deliver therapy, 
manage therapy, or adjust therapy. For example, the historical record of the cardiac 
signal, stored as a sequence of characteristic points, may reveal that a particular 
therapy regimen is appropriate for a particular patient. The long term memory 
capacity of the present subject matter enables storage of data for extended periods of 
time. 

Conclusion 

The above description is intended to be illustrative, and not restrictive. 
Many other embodiments will be apparent to those of skill in the art upon reviewing 
the above description. 
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